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EXPERIMENTAL EVIDENCE OF SUSTAINED COUPLED LONGITUDINAL 
AND LATERAL OSCILLATIONS FROM A ROCKET-PROPELLED MODEL 
OF A 55° SWEPT WING AIRPLANE CONFIGURATION 
By James H. Parks 


SUMMARY 


A rocket-propelled model of a representative 35° sweptback wing 
airplane design has been flight tested at high subsonic Mach numbers to 
investigate the possibility of the existence of coupled longitudinal and 
lateral motions. The model contained pulse rockets arranged to produce 
impulses in yaw. The flight time history is presented herein with some 
semi quantitative analysis. 

The data indicate the existence of sustained longitudinal motions 
of appreciable amplitude forced by lateral oscillations. These longi- 
tudinal motions had twice the frequency of the lateral oscillations and 
are shown to be the result of aerodynamic moments due to sideslip and 
inertial cross coupling. 

Feedback of the longitudinal motions into the lateral character- 
istics through further inertial cross coupling is shown. The roll char- 
acteristics are indicated to be the predominant influence in the inertial 
cross coupling. 


INTRODUCTION 


The concept of a pitching moment arising from a change in sideslip 
angle due to changes in flow conditions at the horizontal tail has been 
recognized for some time, and its effect is noted in reference 1. Other 
investigations (ref. 2, for example) have indicated the possibility of a 
more complex pitching moment arising from inertial coupling of yawing and 
rolling velocities. Recent rocket-model investigations have tended to 
confirm this possibility . 
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To ascertain some idea of the character and nature of these coupled 
longitudinal and lateral motions, the Langley Pilotless Aircraft Research 
Division launched a rocket model of a representative 35° swept wing air- 
plane equipped with pulse rockets arranged to produce disturbances in yaw. 
This preliminary model contained only the minimum instrumentation neces- 
sary to record such motions; consequently, results are presented as 
experimental evidence of sustained coupled longitudinal and lateral 
motions with some observations and semi quantitative analysis. 


SYMBOLS 


M 

V 

W 

s 

b 

c 

g 

a y 

a z 

C N 


C C 

Cy 

C L 

C D 

Cm 


Mach number 
velocity, ft/sec 
weight, lb 
wing area, sq ft 
wing span, ft 

mean aerodynamic chord of wing, ft 
acceleration due to gravity, ft /sec 2 
normal acceleration, ft/sec 2 
lateral acceleration, ft/sec 2 
longitudinal acceleration, ft /sec 2 

normal-force coefficient, 

chord -force coefficient, 

side -force coefficient, — 1 

Sq g 

lift coefficient, Cjj cos a - C c sin a 
drag coefficient, Cq cos a + Cjj sin a 
pitching -moment coefficient 


W_ fn 

Sq g 

. X tk 

Sq g 
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q. 

p 

t 

e 

n 

a 

3 

^max 

9 

9 

7 



yawing -moment coefficient 
dynamic pressure, 0 . 7pM^ , lb/sq ft 
local static pressure, lb/sq ft 
flight time, sec 

inclination of principal longitudinal axis with respect 
to the fuselage center line, deg 

inclination of principal longitudinal axis with respect 
to flight path, deg 

angle of attack, deg or radians 

angle of sideslip, deg or radians 

maximum angle of sideslip for particular oscillations, 
deg 

inclination of X-axis to horizontal, radians 
angle of bank, radians 
angle of yaw, radians 
flight-path angle, radians 

moment of inertia about principal longitudinal axis, 
slug-ft^ 

p 

moment of inertia about principal lateral axis, slug-frt 


moment of inertia about principal vertical axis, 
slug -ft ^ 

moment of inertia about X-axis, I^ cos q + I^ sin t), 
slug-ft^ 

o 

moment of inertia about Y-axis, Iy , slug -ft 

moment of inertia about Z-axis, I 7 cos n + I Y sin ti, 
o -Nd 

slug-ft 

product of inertia, ~(^z “ ^X ) sin T ' cos ^ slug-ft^ 
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X 

Y 

Z 

L 

M 

N 


force along X-axis, lb 
force along Y-axis, lb 
force along Z-axis, lb 
rolling moment, ft -lb 
pitching moment, ft -lb 
yawing moment, ft -lb 


a, p, 9, ijf, cp 


0, cp 


rotational velocity; for 

rotational acceleration; 
radians /sec 2 


example, 0 = — , 
dt 

for example, 0 = 


radians/sec 



a = i(-Si - ^ + ^V -kJ 

2 \m' I Y ' /\ 57 * 3 / 


00 

k 


forcing frequency 

natural frequency of the model. 



Sx C ^V 1 \ 

I Y ' m'ly' /\ 5 T - 3 / 


constant factor used to define a parabolic distribu- 
tion of aerodynamic pitching moment with p 


Cj^, Cj^, C^*, rate of change of aerodynamic forces and moments 

with respect to subscripts; for example, 

^ da 


m 

m' 


R 


mass of the model 

mV 

qS 

3l 

qSc 

Reynolds numbers, based on c 
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MODEL AND APPARATUS 


A three -view drawing of the model is shown in figure 1 and the 
geometric and mass characteristics are tabulated in table I. Photographs 
of the model are shown as figure 2. Details of the construction of this 
type model and the basic test technique are discussed in reference 3- A 
discussion of the pulse-rocket technique is presented in reference L. 
Locations of the pulse rockets are indicated in figure 1 and one exhaust 
port is evident in figure 2(a). 

The model contained a six-channel telemetering installation which 
transmitted continuous records of normal, transverse, and longitudinal 
accelerations, angle of attack, angle of sideslip, and total pressure. 
Tracking radar and radiosonde data were used to obtain local static 
pressure during flight. 

It is believed that the errors with this type instrumentation are 
primarily systematic in nature and the absolute levels are within ± 0.007 
in C^, ±0.008 in Cy, and ±0.003 in Up at M = 0.80. Incremental 
values or slopes and variations with Mach number of the various quanti- 
ties should be more accurate than the absolute values. In general, the 
individual instruments had natural frequencies of the order of 120 cps 
and damping ratios of 0 . 7 . The accuracy of the combined a and (3 
instrument is not completely known; however, similar vane-type instru- 
ments sensitive to angle of attack only have been shown to be accurate 
to about 0.3°. The Mach number is believed correct within 0.02. 

The measurements made to determine the inclination of the principal 
longitudinal axis with respect to the fuselage center line are not con- 
sidered accurate enough to determine values of less than 1°. Therefore, 
the value of e = 0° shown in table I probably is some small negative 
angle of less than 1°. 


RESULTS AND DISCUSSION 


The axis system used in this discussion is shown in figure 3» The 
positive directions of forces, moments, and angles are indicated in the 
figure. Time histories of Mach number, dynamic pressure, and Reynolds 
number are shown in figure L. Angle of attack, angle of sideslip, lift 
coefficient, side-force coefficient, and a drag coefficient are shown 
in figure 5 as a function of flight time. As noted in the figures, the 
normal accelerometer recorded only positive values and, thus, no negative 
lift-coefficient values are available below approximately -0.02. 
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Although no measurement of the rolling velocity of accuracy compa- 
rable to that of the telemetered data was made, special radio equipment 
was employed to provide an order of magnitude of the rolling velocities 
encountered. These data are shown in figure 6. The rolling velocity 
of approximately ±10 radians per second with periods of the order of 
0.5 second indicates that the model reached bank angles of approxi- 
mately ±50°. The order of magnitude of these bank angles was verified 
by the motion-picture history of the flight. 

It should be emphasized here that the phase relationships between 
the lateral motions are quite important in the analysis of the total 
motion and cannot be determined from the data available with sufficient 
accuracy. Thus, a direct analysis of the record is virtually impossible. 

The initial longitudinal disturbance (fig. 5) near t = 3.0 seconds 
is complicated by the firing of the first pulse rocket before the sus- 
tainer rocket motor burned out. At sustainer motor burnout, a longi- 
tudinal oscillation was induced by a change in trim because the thrust 
line was not through the model center of gravity; thus, the a and 
Cl values near t = 3-0 seconds cannot be interpreted as resulting 

solely from coupled motion. 

The Dutch roll motion, as indicated by the sideslip oscillations, 
is poorly damped throughout the flight and appears dynamically unstable 
in the time interval t 83 3*6 to t 83 5*0 seconds. These trends are 
in general agreement with the results of reference 5 wherein dynamic 
lateral instability is indicated for several flight conditions. The 
lateral stability characteristics of similar configurations have been 
dealt with extensively in references 5 to 8 among others and it is beyond 
the scope of this report to analyze the contributing factors to these 
characteristics . 


Aerodynamic Pitching Moment Due to Sideslip 

The frequency relationship shown between a and 3 (period of a 
is one-half the period of 3) considered characteristic of the concept 
of coupled motions arising from changes in flow conditions of the hori- 
zontal tail since a longitudinal response to sideslip should be inde- 
pendent of the direction of sideslip. This effect is clearly shown by 
the wind-tunnel data of figure 7. Obviously, the shape of this curve 
will be somewhat altered at higher Mach numbers and different angles of 
attack; however, for the present discussion, an extrapolation of these 
data is not considered necessary. 
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The longitudinal stability data of reference 9 and the longitudinal 
trim data of reference 3 were used in conjunction with the data of fig- 
ure 7 to approximate the angle -of -attack response by assuming 
3 = 3 max sin cut and C m (p) = kp^. When this forcing function is 
used in the usual two-degrees-of -freedom longitudinal equation in the 


form 


kft mflY 2 5in 2 mt 
a + 2aa 4* a^a = — — j , 

I Y 


it may be shown that the steady-state 


angle -of -attack response will be 


a = a. 


trim 


+ ^Pmax 2 ( 1 


sin 2(cot + 8) 


21-t 




[k 2 - 4cu 2 ) 2 + l6a 2 ou 2 
-1 (% 2 " ik ° 2 ) 


where 


9 = i tan 
2 


172 


4ao> 


Values of a responses computed by this method are compared with 
representative flight values in figure 8 at flight times where the angle- 
of -sideslip variation may be reasonably approximated by a sine wave. The 
forms of the angle -of -attack responses are quite comparable in both the 
computed and actual cases. While the continuous changes in altitude and 
Mach number introduce changes in the aerodynamic forces and moments con- 
sidered as constants in the above equations (see refs. 5 and- 6), these 
aerodynamic variations are not considered powerful enough to produce the 
differences shown which are most noticeable as apparent irregularities in 
the recorded angle of attack. 


Applications of the Equations of Motion 
Including Inertial Cross-Coupling Terms 

It is important to note the rapid increase in a between t » 3*5 
and t *» 4.1 seconds (fig. 5)- The incremental angle of attack is 
approximately doubled in two cycles of the relatively constant amplitude 
yawing oscillation. Also the rate of roll (fig. 6) indicates a decided 
increase at these flight times. Subsequent to the buildup in a, 
nonlinearities appear in the lateral parameters p and Oy. Thus, it 
appears that the nature of the coupling is possibly dependent upon other 
factors in addition to the sideslip. 

Further explanations for these coupled motions may be found by 
examination of the equations of motion for the six degrees of freedom in 
nonlinear form. These equations have been developed by various authors 
(ref. 2, for example) and may be expressed as follows: 
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mV 



+ a9 



= X 


( 1 ) 


mV(p + ijr - a4>) = Y 


( 2 ) 


mV(d + 3<p - 0) = Z ( 5 ) 

IjjT + - 1x2 (t + cpe) = L (4) 

T Y® + ( X X “ X z)^ " ~ 9 2 ) = M ( 5 ) 

V + (h ~ z x)*i ~ x xzfe - ®*) = N (6) 


Note particularly the inertial terms in equations ( 5 ) and ( 6 ) which may 
become appreciable if ly and are large relative to Iy. These 

terms become particularly important where large roll rates are involved 
since cp is a common factor in three terms and cp appears in the 
fourth. In equation (5)* it is indicated that a pitching motion may be 
forced by a pure Dutch roll motion through inertial coupling. 


These conditions are fulfilled in the present tests; that is, 

! y » I z » I x (table I) and high roll rates are indicated in figure 6 . 


Typical values of the basic rotational velocities are shown in 
figure 9 at flight times where these values are near maximums. These 
data were determined as follows: 0 = d + 7 where d is from differ- 

entiation of angle-of -attack— time data and y « £E£L S ., y wa s computed 

by differentiation of angle-of-sideslip — time data and using the assump- 
tion '(r = -|3, and the data of figure 6 were faired to obtain a contin- 
uous 9 . It is readily apparent from the relative magnitude of these 
velocities that the roll rate is the predominant factor and should be 
known at least as accurately as the other data. 


These velocities were used to compute the cross-coupling inertial 
forcing functions which are shown in figure 10. The terms I ^ 2 and 

I XZ 0 ^ are omi tted because 9 » if 2 and 9 » 0 \jr. Included in the 

pitching -moment summation is the aerodynamic moment arising from sideslip 
as discussed previously. 
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For the longitudinal case, it is shown that the inertial moments are 
of approximately the same frequency, but of opposite # sign from the aero- 
dynamic moment due to sideslip, with the (l x - I z )<pt term generally 
having twice the amplitude of the Ixz^ term. The summation of all 
three terms indicates an input generally having the form of p max sin(2ajt) 

similar to that indicated for the aerodynamic term alone, but of con- 
siderably greater magnitude. 

For the lateral, case, the inertial terms generally are of the same 
order of magnitude but with the I^qj term effectively leading the 

(iy - l x )0<p term. The summation of the two effects indicates a yawing 
moment generally having the form of sin(3mt) and of larger ampli- 

tude than either of the contributing effects. 

Harmonic Analyses of the Responses 

The foregoing concepts can be further studied by harmonic analyses 
of the responses. The twelve ordinate scheme of harmonic analysis was 
used to determine the predominant frequencies contained in the angle-of- 
attack and angle -of -side slip data. 

The Fourier series representing the angle -of -attack data was calcu- 
lated to be 

a = -1.71 - 0.30 sin(ajt - 85 °) + 2.40 sin(2ajt + 79.2°) - 
O .38 sin(3ajt - 64.8°) + 0.13 sin(4cot - 18 °) + 

0.08 sin(5cut + 45°) + 0.05 sin(6ajt - 90°) 

This equation is represented graphically in the top plot of fig- 
ure 11. Only the first three harmonics are shown, since by inspection 
they are of primary importance. As indicated previously, those inputs 
having a frequency of 2ajt account almost entirely for the angle of 
attack. It may be of interest to note that as shown by the con- 
stant -I. 7 I, the average effect of the various pitching moments is to 
induce an effective change in trim of about -1.2° since the longitudinal 
trim alone was determined as about - 0 . 5 ° from reference 3 . 

The Fourier series representing the angle of sideslip was calculated 

to be 

0 = 0.12 - 7.57 sin(cot + 58 . 2 °) + 0.88 sin( 2 cut + 52.4°) + 

0.90 sin(3aft - 48.2°) + 0.12 sin(4ojfc) + 0.08 sin(5ojt + 60 . 5 °) - 
0.03 sin(6ctft + 90°) 
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As in the longitudinal case, the first three harmonics are represented 
graphically. (See middle plot of fig. 11.) The angle-of -sideslip oscil- 
lation is composed almost entirely of the fundamental harmonic. The 
shape is distorted by the second and third harmonics which are approxi- 
mately equal in magnitude. It would be deduced from the basic nonlinear 
equations that the inertial yawing moments would have the form of 
sin(3<x>t ) . The presence of an appreciable second harmonic is believed 
indicative of the fact that these inertial terms are not true third 
harmonics of the fundamental frequency. 

The third harmonic is more easily identified in the side-force data. 
The respective Fourier series is 

Cy = 0.012 + O.O67 sin(cot + 71*9°) + 0.00L sin(2oot + 63.5°) + 

0.016 sin(3cot + 75. 9°) + 0.001 sin(4<jot) + 0.003 sin(5<nt + 90°) + 
0.001 sin(6o)t + 90°) 

The first three harmonics of these data are shown in the bottom plot of 
figure 11 where obviously the first and third harmonics are predominant. 
The third harmonic represents the a<p contribution as shown in equa- 
tion (2). 


Response to Total Pitching Moment 

Since the harmonic analysis indicates that the angle -of -attack 
response is primarily a second harmonic of the lateral oscillation, the 
response to a pitching moment having the form suggested in figure 11 and 
the amplitude indicated by figure 10 was computed using the two-degree- 
of -freedom longitudinal equations presented earlier. The results are 
compared with the actual response in figure 12. The agreement is con- 
sidered very good considering the approximations used. 


Examination of Roll Equation 

It is interesting to note that the rolling velocity is probably 
relatively unaffected by inertial cross coupling. The term (l z - ly)0ty 
in equation (4) should be small because I z 555 Iy (table I) and both \(r 
and 0 are much smaller than cp. Also, in the product -of -inertia 
term, (ty + <p9 ) is small in comparison with cp or 9^. 
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General Remarks 

While no direct evidence is available, it is believed worth some 
conjecture concerning the buildup in loads. If the resonant frequency- 
in pitch approaches twice the resonant frequency in sideslip, which was 
true in the present tests, it appears that prohibitive lift loads may 
develop rapidly during lateral maneuvers. It should be emphasized that 
a major portion of these loads occurs as a result of velocities and thus 
large amplitude displacements might be tolerated if the periods of the 
lateral motions are sufficiently long. 


CONCLUDING REMARKS 


From the flight time history of a rocket -propelled model of a repre- 
sentative 35° sweptback wing airplane, it is indicated that coupled 
longitudinal motions were excited and sustained by pure lateral oscilla- 
tions. The resulting longitudinal motions had twice the frequency of the 
lateral oscillations and rapidly developed lift loads of appreciable 
magnitude . 

The longitudinal moments are attributed to two sources, aerodynamic 
moments due to sideslip and inertial cross coupling. The roll character- 
istics are indicated to be the predominating influence in the inertial 
cross-coupling terms. 

At the higher pitching rates, irregularities in the lateral data 
are noted. These changes are shown to be consistant with the inertial 
cross -coupling terms in the yawing -moment and side-force equations. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March Jl, 195^- 
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TABLE I 

MODEL GEOMETRY AND MASS CHARACTERISTICS 


Wing: 

Root airfoil section (normal to 0.30 chord) NACA 63-010 

Tip airfoil section (normal to 0.30 chord) NACA 63^-012 

Total area, sq ft 2.97 

Taper ratio 0.57 

Aspect ratio 3. 57 

Sweepback at 0.30 chord, deg 35 .0 

Incidence at fuselage center line, deg 3*0 

Dihedral, deg -3.0 

Geometric twist, deg 0 


Horizontal tail: 

Airfoil section (normal to 0.30 chord) NACA 63-010 

Total area, sq ft 0.68 

Taper ratio 0.50 

Sweepback, deg 46.0 

Aspect ratio 3*58 

Dihedral, deg 0 

Geometric twist, deg 0 


Vertical tail: 

Airfoil section (normal to 0.30 chord) NACA 63-OIO 

Total area (to fuselage center line), sq ft 0.94 

Taper ratio 0.30 

Aspect ratio (based on double span) 3*57 

Sweepback at 0.30 chord, deg 49 


Weight, lb 

Center-of -gravity location, percent c . 
IXq* slug-ft 2 

Iy , SlUg-ft 2 

0 2 

I z , slug-ft 

€ (referred to fuselage center line), deg 


56.56 

14.2 

0.26 

3.21 

3.86 

0 
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Figure 1.- General arrangement of the model. All dimensions in inches. 
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Figure 2.- Photographs of the model. 
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(b) Model plus booster in launching position. 


L-76551+.1 


Figure 2.- Concluded. 
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Figure 3*- The stability system of axes. Arrows indicate positive 
directions of moments, forces and angles. The origin is at the 
center of gravity. 
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/^//g/yf f/rrxz, sec 

Figure 4.- Time history of altitude, Mach number, dynamic pressure, and 

Reynolds number. 
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Figure 5 .- Time history of the flight. 
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Figure 5.- Continued. 
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Figure 5*- Concluded 
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Figure 6.- Rolling velocities as indicated by special radio equipment. 
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Figure 7.- Static variation of pitching-moment coefficient with angle 
of sideslip from wind-tunnel tests. M = 0.17; C T = 0.36* 

a = k •25°; i t = -1.42°. 
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O .01 .06 J 2 ./e .20 .21 26 .22 .36 40 
//?c <2 rr?<2 c? /a/ //me, sec 

(a) Flight time near k.2 seconds. 



O 04 .OS ,/2 ./6 .20 .24 28 .32 .36 .40 .44 .48 .62 .06 

/sTcreman/a/ t//ne , sec 

(b) Flight time near 10.2 seconds. 

Figure' 8.- Comparisons of actual and computed dynamic responses in angle- 
of-attack to angle-of-sideslip variations. 
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/ r V/p'/7^ frmcz, 5<zc 


Figure 9 *- Rotational velocities at a representative flight time. 
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/^/sgr/yf '/‘/mo., sec 


(a) Pitching moments. 



F//ght h/7?<2, 5<zc 


(b) Yawing moments. 

Figure 10.- Inertial cross-coupling moments at a representative flight time. 
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Figure 11.- Harmonic analyses of typical data. 
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f'/rgh / t/nne.^ s<zc 


Figure 12.- Comparison of computed and actual angle-of-attack responses 

to total cross- coupling moments. 
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